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Abstract — In this paper, the design of a wireless communica- 
tion device relying exclusively on energy harvesting is considered. 
Due to the inability of rechargeable energy sources to charge 
and discharge at the same time, a constraint we term the energy 
half-duplex constraint, two rechargeable energy storage devices 
(ESDs) are assumed so that at any given time, there is always 
one ESD being recharged. The energy harvesting rate is assumed 
to be a random variable that is constant over the time interval 
of interest. A save-then-transmit (ST) protocol is introduced, in 
which a fraction of time p (dubbed the save-ratio) is devoted 
exclusively to energy harvesting, with the remaining fraction 1—p 
used for data transmission. The ratio of the energy obtainable 
from an ESD to the energy harvested is termed the energy storage 
efficiency, r\. We address the practical case of the secondary ESD 
being a battery with 77 < 1, and the main ESD being a super- 
capacitor with r\ = 1. Important properties of the optimal save- 
ratio that minimizes outage probability are derived, from which 
useful design guidelines are drawn. In addition, we compare 
the outage performance of random power supply to that of 
constant power supply over the Rayleigh fading channel. The 
diversity order with random power is shown to be the same as 
that of constant power, but the performance gap can be large. 
Finally, we extend the proposed ST protocol to wireless networks 
with multiple transmitters. It is shown that the system-level 
outage performance is critically dependent on the number of 
transmitters and the optimal save-ratio for single-channel outage 
minimization. 

Index Terms — Energy harvesting, save-then-transmit protocol, 
outage minimization, fading channel, energy half-duplex con- 
straint, energy storage efficiency, TDMA. 



I. Introduction 

The operation of communication networks powered either 
largely or exclusively by renewable sources has become 
increasingly attractive, both due to the increased desire to 
reduce energy consumption in human activities at large, and 
due to necessity brought about by the concept of networking 
heterogeneous devices ranging from medical sensors on/in the 
human body to environment sensors in the wilderness J3], 
H). Sensor nodes are powered by batteries that often cannot 
be replaced because of the inaccessibility of the devices. 
Therefore, once the battery of a sensor node is exhausted, the 
node dies. Thus the potentially maintenance-free and virtually 
perpetual operation offered by energy harvesting, whereby 
energy is extracted from the environment, is appealing. 
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The availability of an inexhaustible but unreliable energy 
source changes a system designer's options considerably, com- 
pared to the conventional cases of an inexhaustible reliable en- 
ergy source (powered by the grid), and an exhaustible reliable 
energy source (powered by batteries). There has been recent 
research on understanding data packet scheduling with an 
energy harvesting transmitter that has a rechargeable battery, 
most of which employed a deterministic energy harvesting 
model. In (3), the transmission time for a given amount of data 
was minimized through power control based on known energy 
arrivals over all time. Structural properties of the optimum 
solution were then used to establish a fast search algorithm. 
This work has been extended to battery limited cases in 
ID, battery imperfections in J5), [ED, and the Gaussian relay 
channel in (7). Energy harvesting with channel fading has 
been investigated in |<X| and [9), wherein a water-filling energy 
allocation solution where the so-called water levels follow a 
staircase function was proved to be optimal. 

In scenarios where multiple energy harvesting wireless 
devices interact with each other, the design needs to adopt 
a system-level approach lfl2l . |[T3l . fl4l . In |[T3l . the medium 
access control (MAC) protocols for single-hop wireless sensor 
networks, operated by energy harvesting capable devices, were 
designed and analyzed. In fl4l . N energy harvesting nodes 
with independent data and energy queues were considered, 
and the queue stability was analyzed under different MAC 
protocols. An information theoretic analysis of energy har- 
vesting communication systems has been provided in lfT31 . 
Ifl6| . In |fT31l , the authors proved that the capacity of the 
AWGN channel with stochastic energy arrivals is equal to 
the capacity with an average power constraint equal to the 
average recharge rate. This work has been extended in [TT61 
to the fading Gaussian channels with perfect/no channel state 
information at the transmitter. 

Due to the theoretical intractability of online power schedul- 
ing under the energy causality constraint (the cumulative 
energy consumed is not allowed to exceed the cumulative 
energy harvested at every point in time), most current research 
is focused on an offline strategy with deterministic channel and 
energy state information, which is not practical and can only 
provide an upper bound on system performance. An earlier 
line of research considers the problem of energy management, 
with only causal energy state information, in communications 
satellites [10], which formulated the problem of maximizing a 
reward that is linear in the energy as a dynamic programming 
problem. In ifTTl . energy management policies which stabilize 
the data queue have been proposed for single-user communi- 
cation under linear energy-rate approximations. 
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Fig. 1. Save-Then-Transmit (ST) Protocol 

In this paper, we focus our study on the design of practical 
circuit model and transmission protocol for energy harvesting 
wireless transmitters. To be more specific, we consider a 
wireless system with one transmitter and one receiver, with 
the transmitter using a save-then-transmit (ST) protocol (see 
Fig. [T]) to deliver Q bits within T seconds, the duration of 
a transmission frame. Because rechargeable energy storage 
devices (ESDs) cannot both charge and discharge simultane- 
ously (the energy half-duplex constraint), an energy harvesting 
transmitter needs two ESDs, which we call the main ESD 
(MESD) and secondary ESD (SESD). The transmitter draws 
power from the MESD for data transmission, over which time 
the SESD is connected to the energy source and charges up. 
At the end of transmission for a frame, the SESD transfers 
its stored energy to the MESD. A fraction p (called the save- 
ratio) of every frame interval is used exclusively for energy 
harvesting by the MESD0 The energy storage efficiency, 
denoted by 77, of each ESD may not be 100 percent, and a 
fixed amount of power P c is assumed to be consumed by the 
transmitter hardware whenever it is powered up. The frame 
interval T is assumed to be small relative to the time constant 
of changes in the ESD charging rate (or energy arrival rate). 
The energy arrival rate is therefore modeled as a random 
variable X in Joules/second, which is assumed to be constant 
over a frame. 

Under the above realistic conditions, we minimize the out- 
age probability (to be defined in the next section) over p, when 
transmitting over a block fading channel with an arbitrary 
fading distribution. In this work, we particularize to the case 
where the MESD is a high-efficiency super-capacitor with 
77=1, and the SESD is a low-efficiency rechargeable battery 
with < ?y < 1. Based on the outage analysis, we compare 
the performance between two system setups: the (new) case 
with random power supply versus the (conventional) case 
with constant power supply, over the Rayleigh fading channel. 
It is shown that energy harvesting, which results in time- 
varying power availability in addition to the randomness of the 
fading channel, may severely degrade the outage performance. 
To be concrete, we further consider exponentially distributed 
random power, and show that although the diversity order with 

'Note that the energy source can be connected only to either the SESD or 
the MESD, but not both. 



exponential power is the same as that with constant power over 
the Rayleigh fading channel, the outage probability curve may 
only display the slope predicted by this diversity analysis at 
substantially higher SNRs. 

Finally, we extend the ST protocol for the single-channel 
case to the general case of wireless network with multiple 
transmitters. We propose a time division multiple access 
(TDMA) based ST (TDMA-ST) protocol to allocate orthogo- 
nal time slots to multiple transmitters that periodically report 
to a fusion center. Specifically, we consider two types of source 
data at transmitters as follows: 

• Independent Data: transmitters send independent data 
packets to the fusion center for independent decoding; 

• Common Data: transmitters send identical data packets 
to the fusion center, where diversity combining is applied 
to decode the common data. 

It is shown that for both cases if the number of transmitters N 
is smaller than the reciprocal of the optimal transmit-ratio (1 — 
p) for the single-channel outage minimization, all transmitters 
can operate at their individual minimum outage probability. 
However, as N goes up and exceeds this threshold, the system- 
level outage performance behaves quite differently for the two 
types of source data. 

The rest of this paper is organized as follows. Section [TT] 
presents the system model. Section [ill] considers finding the 
optimal save-ratio for outage minimization and analyzes its 
various properties. Section |IV] compares the outage perfor- 
mance between fixed power and random power. Section [V] 
introduces the TDMA-ST protocol for the multi-transmitter 
case. Section [VTl shows numerical results. Finally, Section [VlTl 
concludes the paper. 

II. System Model 

A. Definitions and Assumptions 

The block diagram of the system is given in Fig. |2] The 
energy harvested from the environment!! is first stored in either 
the MESD or the SESD at any given time, as indicated by 
switch a, before it is used in data transmission. The MESD 
powers the transmitter directly and usually has high power 
density, good recycle ability and high efficiency, e.g. a super- 
capacitor ifTTl . Since the MESD cannot charge and discharge 
simultaneously, a SESD (e.g. rechargeable battery) stores up 
harvested energy while the transmitter is on, and transfers all 
its stored energy to the MESD once the transmitter is off. We 
assume in the rest of this paper that the SESD is a battery 
with an efficiency 77H where 77 G [0, 1]. This means that a 
fraction 77 of the energy transferred into the SESD during 
charging can be subsequently recovered during discharging. 
The other 1 — 77 fraction of the energy is thus lost, due to e.g. 
battery leakage and/or circuit on/off overhead. The reason of 
choosing a single-throw switch (switch a in Fig. I3 between 
the energy harvesting device (EHD) and ESDs is that splitting 
the harvested energy with a portion going to the SESD, when 
the transmitter does not draw energy from the MESD, is not 

2 Wind, solar, geothermal, etc. 

3 In practice, the battery efficiency can vary from 60% to 99%, depending 
on different recharging technologies 1181 . 
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Fig. 2. Energy Harvesting Circuit Model 



energy efficient due to the SESD's lower efficiency. Note that 
at the current stage of research, the optimal detailed structure 
of an energy harvesting transmitter is not completely known 
and there exist various models in the literature (see e.g. fl8], 
(6), 191). The proposed circuit model, given in Fig. [2] provides 
one possible practical design. 

We assume that Q bits of data are generated and must be 
transmitted within a time slot of duration T seconds (i.e., 
delay constrained). In the proposed ST protocol, the save-ratio 
p is the reserved fraction of time for energy harvesting by 
the MESD within one transmission slot. In other words, data 
delivery only takes place in the last (1 — p)T seconds of each 
time slot, which results in an effective rate of i? e ff = nzrjp 
bits/sec. We also allow for a constant power consumption of P c 
Watts by the transmitter hardware whenever it is powered on. 
The combined influence of p, r\ and P c on outage probability 
is quantified in this work. 

Assume a block-fading frequency-nonselective channel, 
where the channel is constant over the time slot T. Over any 
time slot, the baseband-equivalent channel output is given by 



y = h ■ x + n, 



(1) 



where x is the transmitted signal, y is the received signal, 
and n is i.i.d. circularly symmetric complex Gaussian (CSCG) 
noise with zero mean and variance cr^. 

For any frame, the ST protocol (cf. Fig. [TJ is described as 
follows: 

> During time interval (0, pT], harvested energy accumu- 
lates in the MESD, which corresponds to the situation 
that switches b, c are open and a connects to the MESD 
in Fig. |2j 

■ From time pT to T, the transmitter is powered on for 
transmission with energy from the MESD. Since the 
transmitter has no knowledge of the channel state, we 
assume that all the buffered energy in the MESD is used 
up (best-effort transmission) in each frame. Since the 
MESD cannot charge and discharge at the same time, 
the SESD starts to store up harvested energy while the 
transmitter is on. Referring to Fig. |2] c is closed, b is 
open and a switches to the SESD; 

• At time T, the transmitter completes the transmission and 
powers off. The SESD transfers all its buffered energy to 



the MESD within a negligible charging time, at efficiency 
r\. In other words, b is closed and switches a and c are 
open in Fig. |2] 

B. Outage Probability 

It is clear that the energy harvesting rate X is a non-negative 
random variable with finite support, i.e., < X < Ph < 
00, as the maximum amount of power one can extract from 
any source is finite. Suppose fx{%) and Fx{x) represent its 
probability density function (PDF) and cumulative distribution 
function (CDF), respectively. According to the proposed ST 
protocol, the total buffered energy in the MESD at t = pT 
(the start of data transmission within a transmission slot) is 
given by 

E T = X[p + i 1 (l-p)]T. (2) 

Ej_ = x [ p 

i-p 



Denote P = 



V 



as the average total 



power, which is constant over the entire transmission period, 
and P c as the circuit power (i.e. the power consumed by the 
hardware during data transmission), again assumed constant. 
The mutual information of the channel ([1) conditioned on X 
and the channel gain h is (assuming P > P c ) 

(P-P c )\h\ 2 



R T = log, 1 



= log 2 (1 + (P - P C )T) 
(3) 

where T = with PDF / r (-) and CDF F r ( ). 

For a transmitter with energy harvesting capability and 
working under the ST protocol, the outage event is the union 
of two mutually exclusive events: Circuit Outage and Channel 
Outage. Circuit outage occurs when the MESD has insufficient 
energy stored up at t = pT to even power on the hardware 
for the duration of transmission i.e. Et < -P c (l — P)T or 
equivalent P < P c . Channel outage is defined as the MESD 
having sufficient stored energy but the channel realization does 



not support the required target rate i? c ff 



- Q 



(l-p)T 



bits/s. 



Recalling that X € [0, Pjy], the probabilities of Circuit 
Outage and Channel Outage are therefore: 

pcircuit =pr | p< p c | 

Fx [</>(■)] 'd P B >(/>{■) 
1 otherwise. 

pchannel = ft (l + ( P _ Pc)T) < R eSj p > pj 

O-Reff _ I 

= Pr i r < — —,P > P c 



(4) 



P-Pr 



Ph 



f x (x)F r [g(-)}dx if P H ><£(•) 
otherwise. 



(5) 



where g(p,r),P c ) = x[ _^ +v] _ Pc and <f)(p,r],P c ) = zzr^- 
Since Circuit Outage and Channel Outage are mutually exclu- 
sive, it follows that 



Pn 



-pcircuit 
out 



1 rychannel 

T" r OUt 

Fx [</>(■)] + 
f^ ) fx(x)F r [g(.))dx 
1 



if Ph > </>(■) 
otherwise. 



(6) 
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For convenience, we define 

i-Ph 

Pout(p, V, Pc) = F x [<£(•)] + / fx(x)F T [<?(•)] dx (7) 

where P out (p,r], P c ) < 1 and P H > <j>{-). 

Unlike the conventional definition of outage probability in 
a block fading channel, which is dependent only on the fading 
distribution and a fixed average transmit power constraint, in 
an energy harvesting system with block fading and the ST 
protocol, both transmit power and channel are random, and 
the resulting outage is thus a function of the save-ratio p, the 
battery efficiency 77 and the circuit power P c . 

III. Outage Minimization 

In this section, we design the save-ratio p for the ST 
protocol by solving the optimization problem 

(PI) : min. P out 

0<p<l 

i.e. minimize average outage performance P out in © over 
p, for any given 77 <E [0, 1] and P c £ [0, 00). Denote the 
optimal (minimum) outage probability as P* ut (r],P c ) and the 
optimal save-ratio as p*(rj,P c ). Note that p /• 1 represents 
transmission of a very short burst at the end of each frame, and 
the rest of each frame is reserved for MESD energy harvesting. 
p = is another special case, in which the energy consumed 
in frame i was collected (by the SESD) entirely in frame 
(PI) can always be solved through numerical search, but it is 
challenging to give a closed-form solution for p*(rj,P c ) in 
terms of P c and r\ in general. We will instead analyze how 
p*(r),P c ) varies with P c and 77 and thereby get some insights 
in the rest of this section. 

Proposition 3.1: P ou t(p, "Hi Pc) in © is a non-increasing 
function of battery efficiency 77 and a non-decreasing function 
of circuit power P c for p £ [0, 1). The optimal value of (PI) 
PoutiVi Pc) is strictly decreasing with 77 and strictly increasing 
with P c . 

Proof: Please refer to Appendix [A] ■ 
The intuition of Proposition 13.11 is clear: If 77 grows, the 
energy available to the transmitter can only grow or remain 
the same, whatever the values of p and P c , hence P out must 
be non-increasing with 77; if P c grows, the energy available 
for transmission decreases, leading to higher P ou t. 



A. Ideal System: 77 = 1 and P c = 

Suppose that circuit power is negligible, i.e. all the energy 
is spent on transmission, and the SESD has perfect energy- 
transfer efficiency. The condition Pjj > P c /(jir- + rf) is 
always satisfied, and problem (PI) is simplified to 



(P2) : min. 

o< P <i Jo 



Ph 



fx(x)F T 



(2 (i-p)t 



dx 



where the optimal value of (P2) is denoted as P* ut (l, 0), and 
the optimal save-ratio is denoted as p*(l,0). 



Lemma 3.1: The minimum outage probability when 77 = 1 
and P c = is given by 



fx{x)F T 



2 Q/t _ 1 



dx 



(8) 



and is achieved with the save-ratio p*(l, 0) = 0. 

Proof: Please refer to Appendix [B] ■ 
Lemma im indicates that the optimal strategy for a transmit- 
ter that uses no power to operate its circuitry powered by two 
ESDs with 100 percent efficiency, is to transmit continuously^ 
This is not surprising because the SESD collects energy from 
the environment just as efficiently as the MESD does, and 
so idling the transmitter while the MESD harvests energy 
wastes transmission resources (time) while not reaping any 
gains (energy harvested). However, we will see that this is only 
true when there is no circuit power and the battery efficiency 
is perfect. 

B. Inefficient Battery: 77 < 1 and P c = 

When the SESD energy transfer efficiency 77 < 1 and P c = 
0, (PI) becomes 

(2TT^TT _ 1) 



(P3) 



mm. 

0<p<l 



fx(x)F T 



' p 
■ i-p 



n) 



dx 



where the optimal value of (P3) is denoted as P* ut (r],0), and 
the optimal save-ratio is denoted as p* (77, 0). 

Lemma 3.2: When SESD energy transfer efficiency 77 < 1 
and circuit power P c = 0, the optimal save-ratio p has the 
following properties. 

1) A "phase transition" behavior: 



p*M) = o, 
p*M) >o, 



77 £ 



77 £ 



2^-1 

2T (In 2) fit 



0, 



Q 
2T 

-gr 



(9) 



2T(ln2)f 

2) p* (77, 0) is a non-increasing function of 77, for < 77 < 1. 
Proof: Please refer to Appendix [C] ■ 

According to (0, if the SESD efficiency is above a thresh- 
old, the increased energy available to the transmitter if the 
MESD rather than the SESD collects energy over [0, pT] is 
not sufficient to overcome the extra energy required to transmit 
at the higher rate i? e ff over (pT,T}. The result is that the 
optimal p is 0. On the other hand, if 77 is below that threshold, 
then some amount of time should be spent harvesting energy 
using the higher-efficiency MESD even at the expense of 
losing transmission time. Lemma |3?2l quantifies precisely the 
interplay among 7/, Q, T and p. 

We should note here that even though we consider the 
case of having two ESD's, by setting 7/ = 0, we effectively 
remove the SESD and hence our analysis applies also to 
the single-ESD case. According to (0, if we only have one 
ESD, the optimal save-ratio is p*(0, 0), which is always larger 
than 0. This is intuitively sensible, because with only one 
ESD obeying the energy half-duplex constraint, it would be 
impossible to transmit all the time (p = 0) because that would 
leave no time at all for energy harvesting. 

4 Except for the time needed in each slot to transfer energy from the SESD 
to the MESD, which we assume to be negligible. 
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C. Non-Zero Circuit Power: 77 < 1, P c > 

Non-zero circuit power P c leads to two mutually exclusive 
effects: (i) inability to power on the transmitter for the (1— p)T 
duration of transmission - this is when Ph < </>(•) in ©; 
and (ii) higher outage probability if Ph > </»(■) because some 
power is devoted to running the hardware. 

Since ,j Pc - decreases as p increases, its maximum value 

is — . Therefore, if Ph > — , the transmitter would be able to 
recover enough energy (with non-zero probability) to power on 
the transmitter, i.e. p £ [0, 1). If Ph < — , by condition P H < 

Pc ~ 

P Pc , save-ratio p is required to be larger than — — — . In 
summary, 

. If P c < PhV 

Pout= Pant(p,V,Pc), VpS[0,l) 



If Pc > PhV 



Pout " 



Pout(p,V, P c), 



p< 
p > 



l-»7-h 



(10) 



If P c > t]Ph, referring to ([Tot , we may conclude that 



P*(V,Pc) > 



due to the need to offset circuit power 



consumption. If P c < t]Ph, theoretically, the transmitter is 
able to recover enough energy (with non-zero probability for 
all p £ [0, 1)) to transmit. 

Lemma 3.3: For an energy harvesting transmitter with bat- 
tery efficiency 77 and non-zero circuit power P c , 



n 



Ph " 



2% 



1 



2T(ln2)§ 
Proof: Please refer to Appendix [D] 



p*(r?,Pc)>0. 



(11) 



Intuitively, the smaller the circuit power, the more energy we 
can spend on transmission; the larger the battery efficiency is, 
the more energy we can recover from energy harvesting. Small 
circuit power and high battery efficiency suggest continuous 
transmission (p*(r),P c ) = 0), which is consistent with our 
intuition. According to Lemma 13.31 larger circuit power may 
be compensated by larger ESD efficiency (when the threshold 
for 77 is smaller than 1). A non-zero save-ratio is only desired if 
there exists significant circuit power to be offset or substantial 
ESD inefficiency to be compensated. The threshold depends 
on required transmission rate. 

Remark 3.1: It is worth noticing that if we ignore the 
battery inefficiency or set 77 = 1, Lemma 13.31 could be 
simplified as 



Pc. > 



2*(ln2)£ 



2* 



2$ (In 2)* 



P 



H 



p*(l,P c )>0 (12) 



where only circuit power P c impacts the save-ratio. Since 
the MESD and the SESD are equivalent (77 = 1), harvest- 
ing energy using the MESD is not the reason for delaying 
transmission. Instead, p* > when P c is so large that we 
should transmit over a shorter interval at a higher power, so 
that the actual transmission power minimizes P ou t- Circuit 



power similarly determined the fundamental tradeoff between 
energy efficiency and spectral efficiency (data rate) in |[T9l , 
in which it was shown that with additional circuit power 
making use of all available time for transmission is not the 
best strategy in terms of both energy and spectral efficiency. 
In this paper, outage is minimized through utilizing available 
(random) energy efficiently, wherein circuit power causes a 
similar effect. 

IV. Diversity Analysis 

The outage performance of wireless transmission over fad- 
ing channels at high SNR can be conveniently characterized 
by the so-called diversity order |f20l , which is the high-SNR 
slope of the outage probability determined from a SNR-outage 
plot in the log-log scale. Mathematically, the diversity order 
is defined as 



Um logio(-P OM t) 

7^oc log 10 (7) 



(13) 



where P out is the outage probability and 7 is the average SNR. 

Diversity order under various fading channel conditions 
has been comprehensively analyzed in the literature (see e.g. 
ll20l and references therein). Generally speaking, if the fading 
channel power distribution has an accumulated density near 
zero that can be approximated by a polynomial term, i.e., 
Pr(|/i| 2 < e) rj e fe , where e is an arbitrary small positive 
constant, then the constant k indicates the diversity order of the 
fading channel. For example, in the case of Rayleigh fading 
channel with Pr(|/i| 2 < e) w e, the diversity order is thus 1 
according to (TOT) , 

However, the above diversity analysis is only applicable to 
conventional wireless systems in which the transmitter has 
a constant power supply. Since energy harvesting results in 
random power availability in addition to fading channels, 
the PDF of the receiver SNR due to both random transmit 
power and random channel power may not necessarily be 
polynomially smooth at the origin (as we will show later). 
As a result, the conventional diversity analysis with constant 
transmit power cannot be directly applied. In this section, 
we will investigate the effect of random power on diversity 
analysis, as compared with the conventional constant-power 
case. For clarity, in the rest of this section, we consider the 
ideal system with 77 = 1 and P c = 0, and the Rayleigh fading 
channel with E[T] = E[-££] = f| = A 7 . 

From (0 and ©, the outage probability when 77 = 1 and 
P c = is given by 



Pout - Pr |log 2 (l + PT) < Q p)T 



(14) 



Based on Lemma 13.11 the minimum outage probability is 
achieved with the save-ratio p = 0. Therefore, the outage 
probability is simplified a^fl 



p; ut = p r {pr<c} 



Jo 



fp{p)fr{l)d 1 dp (15) 



5 For convenience, P* ut is used to represent Po Ut (l, 0) in the rest of this 
section. 
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where C = 2t — 1 and the last equality comes from the 
assumption of Rayleigh fading channel so the T is exponential 
distributed. It is worth noting that in this case with r\ = 1 
and p = 0, according to (O, the energy arrival rate X and 
the average total power P are identical and thus can be used 
interchangeably. 

Clearly, the near-zero behavior of P* ut critically depends on 
the PDF of random power fp{p), while intuitively we should 
expect that random power can only degrade the outage per- 
formance. We choose to use the Gamma distribution to model 
the random power P, because the Gamma distribution models 
many positive random variables (RVs) Ell . E2l . The Gamma 
distribution is very general, including exponential, Rayleigh, 
and Chi-Square as special cases; furthermore, the PDF of any 
positive continuous RV can be properly approximated by the 
sum of Gamma PDFs. Supposing that P follows a Gamma 
distribution denoted by P ~ Q({3,A p ), then its PDF is given 
by 



fp(p) = 



p 13 x exp 



Agrcs) 



U(p) 



(16) 



where U(-) is the unit step function, T(-) is the Gamma 
function, and (3 > 0, X p > are given parameters. Referring 
to E3l Lemma 2], which gives the distribution of the product 
of m Gamma RVs, the outage probability in ([T5T l can be 
computed as 



1 g2i ( 



A, 



1 

1,0,0 



(17) 



where G(-) is the Meijer G-function Ell . 

Meijer G-function can in general only be numerically evalu- 
ated and does not give much insights about how random power 
affects the outage performance. Next, we further assume that 
the random power P is exponentially distributed (as a special 
case of Gamma distribution with j3 = 1) to demonstrate the 
effect of random power. 

Lemma 4.1: Suppose that P is exponentially distributed 
with mean A. 



p , the channel is Rayleigh fading with E[T] = 



A 7 , and thus the average received SNR 7 



Ap A 7 



p 2 h . The minimum outage probability P* ut , 
system with rj = 1 and P c = 0, is given by 



under an ideal 



p* 



£ 

k=Q 



1 and P c 

C k+i 



(fc!) 2 (fc + l)7 fc + 1 



1 



C 

In — +2ip(k- 
7 



1 



(18) 
represents 



where ip( x ) is trie digamma function 11241 and In 
the natural logarithm. 

Proof: Please refer to Appendix [E] I 
In the asymptotically high-SNB@ regime, we can approxi 
mate P* ut by taking only the first term of ([T81 as 



P* 



C 
1 



1 - In 



C 

1 



24,(1; 



In 7 

1 



(19) 



6 We assume that high SNR is achieved via decreasing noise power cr^, 
while fixing the average harvested energy. 
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Fig. 3. TDMA based ST (TDMA-ST) 



As observed, P* ut decays as 7 -1 In (7) rather than 7 -1 as in 
the conventional case with constant power, which indicates that 
the PDF of the receiver SNR is no longer polynomially smooth 
near the origin. Hence, the slope of P* ut in the SNR-outage 
plot, or the diversity order, will converge much more slowly to 
7 _1 with SNR than in the constant-power case, suggesting that 
random energy arrival has a significant impact on the diversity 
performance. More specifically, we obtain the diversity order 
in the case of exponentially distributed random power as 



lim 



login 7 + logio ( ln 7) 



1 



(20) 



7^00 l O g 10 7 

which is, in principle, the same as that over the Rayleigh 
fading channel with constant power. We thus conclude that 
diversity order may not be as meaningful a metric of evaluating 
outage performance in the presence of random power, as in 
the conventional case of constant power. 

V. Multiple Transmitters 

In this section, we extend the ST protocol for the single- 
channel case to the more practical case of multiple transmitters 
in a wireless network, and quantify the system-level outage 
performance as a function of the number of transmitters in 
the network. 

A. TDMA-ST 

We consider a wireless network with N transmitters, each 
of which needs to transmit Q bits of data within a time 
frame of duration T seconds to a common fusion center 
(FC). It is assumed that each transmitter is powered by the 
same energy harvesting circuit model as shown in Fig. [2] and 
transmits over the baseband-equivalent channel model given in 
(Q3. We also assume a homogeneous system setup, in which 
the channel gains, energy harvesting rates or additive noises 
for all transmitter-FC links are independent and identically 
distributed (i.i.d). 

In order to allow multiple transmitters to communicate with 
the FC, we propose a TDMA based ST (TDMA-ST) protocol 
as follows (cf. Fig. |3J: 
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Every frame is equally divided into N orthogonal time 
slots with each slot equal to S seconds. 
Assuming perfect time synchronization, each transmitter 
is assigned a different (periodically repeating) time slot 
for transmission, i.e., in each frame, transmitter i is 
allocated the time slot 
Assuming pi = p for al 



(j-l) rp J_rp 

N ± ' N 



1 < i < N. 
i's, each transmitter implements 
the ST protocol with the transmission time in each frame 
aligned to be within its assigned time slot; as a result, 
the maximum transmit- ratio, denoted by 1 — p, for each 
transmitter cannot exceed 1/N, which means that p > 

The protocol described above is illustrated in Fig. [3] Unlike 
the single-channel case where the transmitter can select any 
save-ratio p in the interval < p < 1, in the case of TDMA- 
ST, p is further constrained by p > 1 — -h to ensure orthogonal 
transmissions by all transmitters. Due to this limitation, each 
transmitter may not be able to work at its individual minimum 
outage probability unless the corresponding optimal save-ratio 
p* satisfies p* > 1 — 4 or N < rz-r- In this case, ST 



protocol naturally extends to TDMA-ST with every transmitter 
operating at the optimal save-ratio p* . However, if N exceeds 
the threshold 1 „ , transmitters have to deviate from p* 
to maintain orthogonal transmissions. Next, we evaluate the 
performance of TDMA-ST for two types of source data at 
transmitters: Independent Data and Common Data. 

B. Independent Data 

First, consider the case where all transmitters send indepen- 
dent data packets to the FC in each frame, which are decoded 
separately at the FC. Under the symmetric setup, for a given 
p, all transmitters should have the same average outage per- 
formance. Consequently, the system-level outage performance 
in the case of independent data can be equivalently measured 
by that of the individual transmitter, i.e., 



Pout{p,V, P c)- 



(21) 



We can further investigate the following two cases: 
. N< 



i-p* 



In this case, the additional constraint due to TDMA, p* > 
1 — jj, is satisfied. Since P^ ut is the same as that of 
the single-transmitter case, the system is optimized when 
all transmitters work at their individual minimum outage 
with save-ratio p* . Thus, the minimum system outage 
probability is P D s * t = P* ut {rf, P c ). 
. N> T ± F 

In this case, the TDMA constraint on p* is violated and 
thus we are not able to allocate all transmitters the save- 
ratio p*, which means that each transmitter has to deviate 
from its minimum outage point. Since in this case p* < 
1 — i < p, the best strategy for each transmitter is to 
choose p=l-j r . Thus, P Q s * t = P out (l - i,r?, P c ). 
To summarize, the optimal strategy for each transmitter in 
the case of independent data is given by 

l _ 

(22) 



P 



P, 
1 



TV < 



1-p* 



which implies that the number of transmitters should be kept 
below the reciprocal of the single-channel optimal transmit- 
ratio; otherwise, the system outage performance will degrade. 

C. Common Data 

Next, consider the case where all transmitters send identical 
data packets in each frame to the FC, which applies diversity 
combining techniques to decode the common data. For sim- 
plicity, we consider selection combining (SC) at the receiver, 
but similar results can be obtained for other diversity combin- 
ing techniques |20l . With SC, the system outage probability 
is given by 



P s 



(p out ( P ,n,p c )) 



N 



(23) 



Similarly to the case of independent data, we can get exactly 
the same result for the optimal transmit strategy given in (l22l ) 
for the common-data case, with which the minimum system 
outage probability is obtained as 



P* 



{Pout(l-jf,V,Pc)) 



N 



N < 
N > 



i-p* 
l 

i-p* 



(24) 



From the above, it is evident that the system outage probability 
initially drops as N increases, provided that N < j-z^r- 
However, when N > ^ , it is not immediately clear 
whether the system outage increases or decreases with N, 
since increasing N improves the SC diversity, but also makes 
each transmitter deviate even further from its minimum outage 
save-ratio according to 



VI. Numerical Examples 

In this section, we provide numerical examples to validate 
our claims. We assume that the energy harvesting rate X 
follows a uniform distribution (unless specified otherwise) 
within [0, 100] (i.e., P H = 100 J/s), and the channel is 
Rayleigh fading with exponentially distributed T with param- 
eter A = 0.02. We also assume the target transmission rate 
R Teq = 9. = 2 bits/sQ 

Fig. [4] demonstrates how battery efficiency rj and circuit 
power P c affect the optimal save-ratio p* for the single- 
channel case. As observed, larger P c and smaller rj result in 
larger p*, i.e. shorter transmission time. Since the increment 
is more substantial along P c axis, circuit power has a larger 
influence on the optimal save-ratio compared with battery 
efficiency. p*(l,0) = verifies the result of Lemma 13.11 
for an ideal system, while p*(i],0) along the line P c = 
demonstrates the "phase transition" behavior stated in Lemma 
13.21 The transition point is observed to be ?] = 0.541, which 
can also be computed from 

Fig. [5] shows the optimal (minimum) outage probability 
P* ut (r],P c ) corresponding to p* in Fig. |4] Consistent with 
Proposition 13.11 P* ut (r),P c ) is observed to be monotonically 
decreasing with battery efficiency i] and monotonically in- 
creasing with circuit power P c . Again, P c affects outage 
performance more significantly than r\. From Fig. [5] we see 

7 This is normalized to a bandwidth of 1 Hz, i.e. -R req is the spectral 
efficiency in bis/s/Hz. 
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that for a reasonable outage probability e.g. below 0.05, P c 
has to be small and rj has to be close to 1. Our results can 
thus be used to find the feasible region in the rj — P c plane 
for a given allowable P ou t- 

Figs. [6] and [7] compare the outage performance with versus 
without save-ratio optimization. In Fig.|6]we fix the normalized 
circuit power ■p s - = 0.5, while in Fig. [7] we fix the battery 
efficiency ?/ = 0.8. We observe that optimizing the save-ratio 
can significantly improve the outage performance. It is worth 
noting that P ou t has an approximately linear relationship with 
the normalized circuit power pg- as observed in Fig. [7] which 
indicates that P c considerably affects the outage performance 
as stated previously. 

In Fig. [8] the outage probability for an ideal system (rj = 1, 
P c = 0) is shown by numerically evaluating dTTb . By fixing 
the mean value of P as E[P] = 50 J/s and varying /3 
for the Gamma distributed power from 1 to 5, the resulting 



Fig. 7. Outage performance comparison: rj = 0.8 



outage performance is compared with the case of constant 
power. As observed, the outage probability increases due to 
the existence of power randomness. As /3 increases, the outage 
curve approaches the case of constant power. In Fig. |9] we 
also plot the outage probability for the ideal system with 
exponentially distributed power based on the approximation 
given in as well as for a non-ideal system with the 

normalized circuit power g^pr =0.1 and battery efficiency 
rj = 0.8. In comparison with the constant-power case, for the 
case of ideal system we observe that the high-SNR slope or 
diversity order with random power clearly converges much 
slower with SNR, which is in accordance with our analysis in 
Section [TV] Furthermore, at P out = 10~ 3 , there is about 10 
dB power penalty observed due to exponential random power, 
even with the same diversity order as the constant-power case. 





1 o~ 4 : -H- Ideal system, constant power 

Ideal system, exponential distributed power 
-0- Ideal system, approximation to exponential power 
-0- Non-idea system, exponential distributed power 



-10 10 20 30 40 50 

7 (dB) 

Fig. 9. Outage probability comparison for ideal (ij = 1, P c = 0) versus 
non-ideal (?? = 0.8, P c = 0.1 * E[P]) systems 



It is also observed that there is a small rising part for the outage 
approximation given in (1191 , since this approximation is only 
valid for sufficiently high SNR values (7 > 10 dB). Finally, 
it is worth noting that the outage probability for the non-ideal 
system eventually saturates with SNR (regardless of how small 
the noise power is or how large the SNR is), which indicates 
that the diversity order is zero for any non-ideal system. 

Fig. [TO] shows the outage performance for the case of 
multiple transmitters operating under the TDMA-ST protocol. 
We set the normalized circuit power -p^ = 0.5 and the 
battery efficiency 77 = 0.9. Then, the optimal save-ratio p* 
for single-transmitter outage minimization can be obtained as 
0.7930 by numerical search. Therefore, the threshold value 
for N in the optimal rule of assigning save-ratio values in 



d22b is 1 } p , = 4.83. For the case of independent data, it is 
observed that when N < 4, the system outage probability 
is constantly equal to the optimal single-transmitter outage 
probability P* ut (0.9, 0.5P#); however, as N > 4, the outage 
probability increases dramatically. In contrast, for the case of 
common data, it is observed that the system outage probability 
decreases initially as N increases, even after the threshold 
value and until N = 7, beyond which it starts increasing. 
This implies that there is an optimal decision on the number 
of transmitters to achieve the optimal outage performance. 

VII. Conclusion 

In this paper, we studied a wireless system under practical 
energy harvesting conditions. Assuming a general model with 
non-ideal energy storage efficiency and transmit circuit power, 
we proposed a Save-then-Transmit (ST) protocol to opti- 
mize the system outage performance via finding the optimal 
save-ratio. We characterized how the optimal save-ratio and 
the minimum outage probability vary with practical system 
parameters. We compared the outage performance between 
random power and constant power under the assumption of 
Rayleigh fading channel. It is shown that random power 
considerably degrades the outage performance. Furthermore, 
we presented a TDMA-ST protocol for wireless networks 
with multiple transmitters. In particular, two types of source 
data are examined: independent data and common data. It is 
shown that if the number of transmitters is smaller than the 
reciprocal of the optimal transmit-ratio for single-transmitter 
outage minimization, each transmitter should work with its 
minimum outage save-ratio; however, when the number of 
transmitters exceeds this threshold, each transmitter has to 
deviate from its individual optimal operating point. 

There are important problems that remain unaddressed in 
this paper and are worth investigating in the future. For 
example, we may consider the effect of different config- 
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urations of battery/supercapacitor and MESD/SESD on the 
system performance. It is also interesting to investigate the 
information-theoretic limits for the ST protocol in the case of 
multiple transmitters using more sophisticated multiple- access 
techniques other than the simple TDMA. 

Appendix A 
Proof of PropositionI3.1I 



According to the Fundamental Theorem of Calculus 
we can derive the first derivative of P ou t(p, rj, P c ) in Q with 
respect to 77, P c and p as 



dP„. 



P- 



Oil 



fx 



p 



P 



p 
l-p 



I, 



P 



p 
i-p 



n 



F r (oo) 



dp 



8/3 



fx(x)fr[g(-)] 



dx 



(25) 



2 (1-P)T _ X 



where /3 could be 77, P c or p and g(p, 77, P c ) = ^ p< 

9g{p,ri,Po) 



i-p 

< (XV77 e [0,1] 



It is easy to verify that - 

and 99(p ^ P f c) > 0,yP c € [0,oo]. Therefore P out (p,V,Pc) 
is strictly decreasing with battery efficiency 77 and strictly 
increasing with circuit power P c . Next, we are going to prove 
the monotonicity of P out and P* ut with battery efficiency 77, 
where circuit power P c is treated as constant. 

The condition Ph > P Pc - in (|6]l could be expressed in 

terms of battery efficiency: 77 > p^ — j^—, then 
p 

L out ' 



1. 

Pout 1 



„ < -Sl _ _e_ 

'/ ^ Pff i_p 

'/ > Ph i_p 



(26) 



Consider the following two cases: 

• Suppose 771 < 772 and p^ — j^— > 

- K j% ~ < m < rn, then 

Pout(p,Vi, p c) = ^Pout{p,m,Pc) and 
Pout(p,V2,Pc) = Pout(p,V2,P c )- Since 
Pout(p,ViPc) is strictly decreasing with battery 
efficiency 77, we have 

Pout{p,m,Pc) > Pout{p,V2,Pc)- 

- If 7ft < j£ - < r? 2 , then P out (p,r]i,Pc) = 1 
and P ut{p,m,Pc) = Pout{p,m,Pc)- Therefore 

Pout(p,Vi,Pc) = 1 > Pout(p,m,Pc)- 

- If 771 < 772 < |jj - j^, then P out (p,Vi, Pc) = 
Pout(p, V2,Pc) = 1. which means 

-Pout (yO, 7?1 j -Pc) = -Pout (/», 7/2 , Pc) ■ 





Pi 


Pc 


P2 


Ph 


1-pf. 


Ph 


l-p* 


Pc 


P* 


Pc 


P*2 


Ph 


1-pl' 


Ph 


l-p* 



• Suppose 771 < 772 and p^ — -^— p < 0, we have p^ — 
< 771 < 7/2- Then it could be easily verified that 

Pout(p,Vl->Pc) > Pout(P,V2,Pc)- 

Combining all the above cases, we can conclude that 
Pout(p, f]-, Pc) is a non-increasing function of battery efficiency 
77 given any non-zero circuit power P c for p £ [0, 1). Next, 
we proceed to prove the monotonicity of P* ut (r], P c ). 

Assuming 771 < 772 again, then we could argue that 
-5^ — < ?7i and -§= -r^-r < 779, where p% and 

Ph 1— p* n Ph I-P2 1 

P2 are the optimal save-ratio for 77 = 771 and 77 = 772, 
respectively. Therefore we only need to consider two cases: 
max -I -Ift- - -|ft - -r%r [■ < 771 < 772 and 771 < 

< 772. From the arguments 
we have given for the proof of the monotonicity of P out we 
know that, under these two conditions we have 

Pout(P,Vl,Pc) > Pout(P,V2,Pc)- 

Therefore, 

PLtiVuPc) > Pout(p*l,r)2,Pc) > P*out(m,Pc) 

which completes the proof of the monotonicity for Po Ut (ri, P c ). 
With similar arguments, we could get the results regarding 
circuit power P c . Proposition 13. II is thus proved. 

Appendix B 
Proof of Lemma[3J] 

Since Fr(-) is non-negative and non-decreasing, we have 
a < b => F r Q < F r f ~J (27) 
for any x £ [0,Pfr]- Since fx(') is non-negative, this leads to 

a<b => J fx(x)Fr(^jdx 

< J " fx{x)Fv {^jdx. 

Given the form of P ou t in Problem (P2), with p appear- 
ing only in the numerator of the argument of Fr(")> we 
conclude that P out is a non-decreasing function of g(p) = 
^2 n-p)T — lj (1 — p). Hence minimizing g(p) is equivalent 
to minimizing P ou t- The first and second derivatives of g(p) 
are 



g'{p) = 2< 1 -f> T (In 2) 
g"{p) = 2T^W(\n2) 



(l-p)T 

2 Q 2 



2<i-p)r + 1 



> since Q > 0. 



T 2 (l - pf 

Let h(p) = g'(p)- From the second equation above, h(p) is 
an increasing function. In the range < p < 1, h(p) is thus 
minimized at p = 0, i.e. the minimum of g'(p) is h(0), given 
by 

g' min = 2T(ln2)|-2* + l (28) 
= 2t fln2 Q/T - l) + 1 > (29) 
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for Q > 0. In other words, the smallest value that the gradient 
of g(p) can take in the range < p < 1 for any feasible Q is 
positive, which implies that g(p) is increasing and therefore 
minimized at p = 0, as claimed. The proof of Lemma 13.11 is 
thus completed. 



Appendix C 
Proof of Lemma[3721 



To prove Property 1, we observe that as noted in the proof 
of Lemma 13.11 P out is a monotonic function of g(p) = 



~{ ^ I rjT ^~ m Problem (P3), hence minimizing g(p) leads 
to the same solution as minimizing P ou t- The first derivative 
of g{p) is 



P^&r (In 2) TT ^ T [p + 77(1 -p)]-2 <^> T + 



9'(p) 



1 



[p + r]{l-p)f 



u(p) 



[p + r)(l-p)}< 



It is clear in the above that the sign of g'(p) is the same as 
that of u(p). Since u(l) = +oo and u(p) is a differentiable 
function, if u(0) is negative then there exists a value p c G 
(0, 1) such that u(p c ) = = g'{p c )- It is easily verified that 
u'(p) > 0; hence p c is the unique optimal value of p in this 
case. Conversely, if there exists an p c such that u(p c ) = 0, 
then u(0) must be negative. Hence u(0) < is a necessary 
and sufficient condition for the optimal p to lie in (0, 1). 

The condition w(0) < translates into the following 
condition on 77, which proves the first part of the lemma: 



t(0)<0 => 2t(1h2)^-2t 



1 < 



7] < 



Q 

2t 



1 



2T(ln2)f 



(30) 



where /3 could be 77, P c or p and g(p,r],P c ) 
Furthermore, we have 



Q 

2 {1-P)T _l 



dg(p) 2T^r(In2) Tr a jT [p + ^(l-p)-(l-p)^ 



dp 



[p + r)(l-p)-(l- P y 



Q 

2(i-p> T - 1 



x[p + ri(l-p)-(l-p)%] 

Kp) 



With similar arguments about u(p) in the proof of Lemma 



13.21 we claim that v(0) < 0,V.t S ( (i Pc | - , Pg] is a sufficient 

condition of having p*(j],P c ) > while P c < r\Pu- 

Since w(0) is an increasing function of x, the condition 

v(0) < 0,Vx G ( /°, , Pgl translates into the following 

i-p ' r i 
condition on r] and P c 



«(0) = 2*(ln2)|^-^)-2*-F 



<2*(In4(^-g)-2* 



1 < 



P c 

< 77- ^ < 



2?-l 



ftr 2T(ln2)g 

— -'7 

Combined with the fact that p*(?7, P c ) > Fff P when P c > 
77PH, we may conclude 



p*(77,P c )>0, 77- ^< 



2#-l 



^ 2T(ln2)^' 



(31) 



Lemma 1331 is thus proved. 



To prove the second point, suppose pl(r]i,0) and /^(^jO) 
are optimal save-ratios of (P3) for SESD efficiencies 771 and 
772, where 771 < 772. Then, u(pl,rji) = and u(p2,i]2) = 0. 
Since 771 < 772 and u(p,rj) is an increasing function of 77, we 
have u(pl,r]2) > 0. Combining what we have that u(p,r]) is 
an increasing function of p, u(p2,r]2) =0 and u(p*, 772) > 0, 
we may conclude /^(^O) < p*{r)\,Q). Lemma [3~2l is thus 
proved. 



Appendix D 
Proof of Lemma[3T3] 



Appendix E 
Proof of Lemma|4~T1 



Let Z = Pr, where P and V are exponential random 
variables with mean A p and A 7 respectively. Then the PDF 
of Z could be derived as follows, 

F z (z) = Pi {PT<z} 



1 / e pA t e Ap dp 

\ Jo 



(32) 



According to the proof of Proposition 13. 11 the first derivative 
of Pout(P) 77, P c ) with respect to 77, P c and p is, 



where Ki(x) is the first-order modified Bessel function of the 
second kind and the last equality is given by l25l §3.324.1]: 



dP ; 



dfi 



fx{x)fv[g{-)]^dx 



9/3 



cxp 
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where M(j3) > 0,3?( 7 ) > 0. Let M = -J — 
derivative of F(z) yields 

f(z) = M j-^^i (2MV5) - 2^ (/u (2MVS))' 
= M I — ^Ai (2MVz) 



Taking the Since D 



2M 2 A (2My/z) 
A ( 2 



A p A 7 



(33) 



where d-^l 2 ) 

9z 



-AVi(*) - f A„(z). 
Next, we characterize the outage probability using 
According to dl5l l. we have 



P* out =Vx[PT <C] 



o A p A 7 



A 2 



Ap A-y 



(34) 



Let X 



and D 



We then have 



,13 

P; ut = 2 K (2Vx) dx. (35) 
Jo 



Using the series presentation 



38.447.3], we have 



K (x) 



.,2/.' 



fc=0 



2 2fe (fc!) 2 



V^fc + 1) (36) 



with the series expansion for the modified Bessel function 
given by 



,.2A- 



2 2fc (fc!) 2 



(37) 



(l35T l could be expanded as 



p* 



where 



E 

k=0 



(kiy 



x k \nxdx + ^(k + l) / x k dx 



D 



is the digamma function 
could be evaluated as 
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— lnlYx) = F ^ 
da: T(x) 

Since the two integrals in 



(38) 
(39) 



~ k dx = 



D k+i 



X 



k In xdx = x 



fc+i 







In . 



1 



= D k+i 



k + 1 (fc + 1) 2 
InD 1 



x=D 



x=0 



k + 1 (fc + l) 2 
where lim x -)-o(a;liia;) = 0. Then we have 



E 

fe=0 



2 

(fc!) 2 fc + 1 



2 V fc + l 



(40) 



thus proved. 



~, ([18]) follows. Lemma |4~T| is 
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